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In contrast to the well-known destabilization of globular proteins by high pressure, re- 
eent work has shown that pressure stabilizes the formation of isolated a-heliees. However 
all simulations to date have obtained a qualitatively opposite result within the experimen¬ 
tal pressure range. We show that using a protein foree field (AmberOSw) parametrized 
in eonjunetion with an aeeurate water model (TIP4P/2005) reeovers the eorreet pressure- 
dependenee and an overall stability diagram for helix formation similar to that from experi¬ 
ment; on the other hand, we eonfirm that using TIP3P water results in a very weak pressure 
destabilization of heliees. By earefully analyzing the eontributing faetors, we show that 
this is not merely a eonsequenee of different peptide eonformations sampled using TIP3P 
Rather, there is a eritieal role for the solvent itself in determining the dependenee of total 
system volume (peptide and solvent) on helix eontent. Helieal peptide struetures exelude a 
smaller volume to water, relative to non-helieal struetures with both the water models, but 
the total system volume for helieal eonformations is higher than non-helieal eonformations 
with TIP3P water at low to intermediate pressures, in eontrast to TIP4P/2005 water. Our 
results further emphasize the importanee of using an aeeurate water model to study protein 
folding under eonditions away from standard temperature and pressure. 
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I. INTRODUCTION 


The dependence of protein folding equilibria on thermodynamic control variables can provide 
important insights into the fundamental forces which stabilize folded structures'^. Variation of 
temperature is the most commonly used approach due to the ease with which this can be achieved. 
A detailed comparison of the temperature-dependence of protein stability with that of the hy¬ 
drophobic effect, for example, strongly supports the role of hydrophobic interactions in stabilizing 
protein structure^^’^*. 

A variable which is less exploited is pressure, due to the greater demands of the experiments 
required^^’"^°. However, it is well-known that high pressure tends to destabilize protein native 
structures"^^’"^^, indicative of a positive change of reaction volume for folding. A possible resolution 
of this initially counterintuitive effecU^ was proposed to be the existence of cavities within the 
folded protein"^"^’"^^, with substantial support for this hypothesis coming from experiments on cavity¬ 
forming mutants"^^’"^^. It was also proposed that high pressures can lead to water penetration of 
protein’s hydrophobic core due to reduced solvent-solute interfacial free energy"*^. 

However, the pressure-dependence of the folding equilibria for small, independently folding 
elements of secondary structure such as hairpins and helices clearly cannot fit the same picture. 
Since they contain no evident internal cavities or well-defined hydrophobic core, any pressure 
dependence would have to come from other effects, which might be obscured when studying the 
overall folding of a globular protein. Recent experimental work on the pressure dependence of the 
helix-coil equilibrium, using either FTIR spectroscopy"^^ "^^ or triplet state quenching experiments^® 
has in fact found the opposite trend to that for protein folding: namely a negative reaction volume 
for helix formation, resulting in pressure stabilization of helices for all positive pressures. The 
origin of this effect is not clear, especially considering the small magnitude of the volume change 
(~ 0.2-1.0 cm^.mol”^ per residue). Molecular simulation could potentially help to explain the 
origin of this result; however, simulation studies of the pressure dependence of helix formation 
have qualitatively contradicted experimental results, finding instead helix destabilization at low to 
intermediate pressures, only turning over to stabilization at very high pressure^ 

Here, we investigate the pressure dependence of helix formation for a model 15-residue helix¬ 
forming peptide using two different force field combinations: the Amber ff03* protein force 
field^"^ together with explicit TIP3P water^^ and the Amber ff03w protein force field^® with the 
TIP4P/2005 water model^^. In agreement with earlier studies with TIP3P water on the effect of 
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pressure on helix formation, we find a positive reaction volume for helix formation at low to in¬ 
termediate pressures. In contrast, a qualitatively correct result is obtained with TIP4P/2005, i.e., 
a negative reaction volume for helix formation. This difference, obtained for the same sequence 
with almost identical protein force fields, suggests a key role for water. There are essentially two 
ways in which water can be envisaged to influence the reaction volume: (i) it may alter the confor¬ 
mations sampled for a given total number of helical residues, particularly for the non-helical states 
and (ii) different water models may be more or less closely coordinated with helical than with 
non-helical structures. While there is no doubt that water influences the conformational sampling, 
we show that in fact the different solvation of helical and non-helical states plays a key role in 
determining the volume changes, and for both water models opposes the much larger decrease in 
volume excluded to water associated with helix formation. The result obtained with TIP3P water 
at low to intermediate pressures is qualitatively inconsistent with experiment, as the total system 
volume is higher for helical structures than for non-helical structures. Our results highlight the 
importance of using an accurate water model for capturing biomolecular equilibria at state points 
away from standard conditions. 


II. METHODS 

A. Simulation methods 

Replica exchange molecular dynamics (REMD) simulations with gromacs 4.0.7 or 4.5.3^^ were 
used to sample the folding of the blocked peptide Ac-(AAQAA) 3 -NH 2 with either (i) the Am¬ 
ber ff03* force field^"* and TIP3P water model^^ or (ii) the Amber ff03w force field^^ and the 
TIP4P/2005 water model^’. For ff03*, 42 replicas spanning a temperature range from 275 to 496 
K were used, and for ff03w, 40 replicas spanning 250 to 452 K. Periodic boundary conditions 
with a truncated octadron minimum image cell of initial edge 4.5 nm were used. Simulations were 
maintained at a constant pressure with a Parrinello-Rahman barostat^^ and temperature was con¬ 
trolled via Langevin dynamics with a friction coefficient of 0.2 ps“^ (ff03*) or 1.0 ps“^ (ff03w). 
Note that the difference in friction coefficients will not alter the equilibrium properties considered 
here, and both values are sufficient to maintain constant temperature. Simulations were run for 
100-200 ns per replica. A seperate set of REMD runs was performed for each pressure, using 1 
bar, 2 kbar, 4 kbar, 8 kbar and 12 kbar for Amber ff03* and 1 bar, 1 kbar, 4 kbar, 8 kbar for Amber 
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ff03w. 


To determine the change of volume associated with helix formation, initial configurations were 
obtained by random selection from the 298 K replica of the ff03* simulations. For each number of 
helical residues (0, 3,4,... 15), 20 configurations were chosen with that helix content. From each 
of these configurations, 10 ns simulations were performed using position restraints to keep the 
peptide close to its initial structure, with force constants of 1000 kJ.mol”^nm“^ on each cartesian 
coordinate, and using both the ff03* and ff03w force fields. Average total system volumes were 
computed from each of these runs. 


B. Helix formation 

We define helical states using the backbone Ramachandran angles, in the spirit of the Lifson- 
Roig theory. We define the helical region of the Ramachandran map as 0 G [—100°, —30°] and 
Ip G [—67°, —7°]. A helical segment is defined as three consecutive residues with their backbone 
torsion angles lying in this range. Therefore the total number of helical residues in the blocked 
15-residue peptide considered can take on values in {0, 3,4, 5,..., 14,15}. This method gives 
results which are consistent with other definitions of helix^^. 


C. Helix-coil models 

Simulation data on helix formation were initially described by a Lifson-Roig model^*’, with nu- 
cleation parameter v and elongation parameter w. These parameters were fitted by maximizing the 
likelihood of the observed conformations, given the equilibrium probabilities generated from the 
Lifson-Roig partition function, using a previously described procedure^"^. Since the experimental 
data had been fitted to the Zimm-Bragg model for helix formation®^ we also converted the fit 
parameters to the Zimm-Bragg nucleation and elongation parameters a and s respectively, using 
the approximate relations^^: 


a 


{i + vY^ 

w 


1 -I- V 


( 1 ) 
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D. Thermodynamic Model 


The dependence of the helix elongation free energy, AGei on pressure and temperature was 
fitted to a thermodynamic model: 


AGei(P,T) = AHo + AC'p(T - Tq) 

-TASo-TACpln{T/To) 

+AV{P - Po) + - Pof 

+Aa{T-To){P-Po). (2) 

In this expression, AHq, AS'o, and AVq are the change of enthalpy, entropy and volume as a result 
of helix elongation at reference conditions, taken to be Tq = 298 K and Pq = 1 bar pressure. In 
addition, a constant change of heat capacity, ACp, change of linear expansion coefficient Aa and 
change of compressibility, Af3 were assumed. The data were fitted by non-linear least squares, 
and errors were estimated by Monte Carlo bootstrapping. 

III. RESULTS 

In order to determine the effects controlling helix stability under pressure, we carefully selected 
two protein force fields. Amber ff03* and Amber ff03w. These models are almost identical, being 
originally based on the standard Amber ff03 force field^^. They only differ in that an additional 
empirically determined Fourier term has been added to the tjj backbone torsion angle in each case, 
in order to approximately match experimental helix propensities near 300 In the case of 

ff03*, the calibration was done in conjunction with the TIP3P water model, while for ff03w it 
was done with TIP4P/2005 water (below, it will be assumed when discussing Amber ff03* and 
ff03w that the TIP3P and TIP4P/2005 water models were used, respectively). This allows us to 
test specifically the effects of the water model using closely related protein force fields that have 
very similar helical populations under standard conditions of temperature and pressure. We study 
the 15-residue peptide Ac-(AAQAA) 3 -NH 2 , as a model system which is known to form helix at 
low temperature^"^, and which has been extensively characterized in previous simulations^"^’^®’^^. 

Replica exchange molecular dynamics (REMD) simulations were performed in order to sample 
the temperature-dependent helix-coil equilibrium. The average fraction helix at each temperature 
is shown in Figure 1, as determined using backbone dihedral angles. Very similar results are 
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obtained using the standard DSSP^^ algorithm, which instead uses backbone hydrogen bonds (See 
Fig. As expected both Amber ff03* and ff03w populate 20 — 30 % helix at 300 K and 1 

bar pressure, where they were parametrized against experimentally determined helix populations. 
For Amber ff03w, all pressures used resulted in a stabilization of the peptide, even at 1 kbar. 
On the other hand, over a wide range, from 1 bar to 8 kbar, pressure had very little effect on the 
overall helix propensity for Amber ff03*: a slight decrease in helix fraction at temperatures greater 
than ~ 300iF is observed. Only for a pressure of 12 kbar is a significant stabilization obtained. 
Qualitatively, these results suggest a negative reaction volume for helix formation using Amber 
ff03w, but a very small or positive reaction volume at low to intermediate pressure using Amber 
ff03*. 

We quantify the effect of pressure on the helix-coil equilibrium using a thermodynamic model. 
However, helix-formation is not a simple two-state process, and involves a broad spectrum of pop¬ 
ulated intermediates. Therefore, we first fit a simple Ising-like statistical mechanics model which 
can capture the helix-coil transition. The two classic partition functions for helix-coil formation 
are those by Zimm and Bragg^^ and by Lifson and Roig^°, which are approximately equivalents^. 
We have determined parameters for both models here. The data were initially fitted to the Lifson- 
Roig model, using a previously described maximum likelihood method^"^, yielding a nucleation 
parameter v and an elongation parameter w. These parameters were then converted into the cor¬ 
responding parameters for the Zimm-Bragg model, a and s. Overall, all these parameters show 
essentially similar trends to the global fraction of helix, but can be more justifiably fitted to a 
thermodynamic two-state model since they describe the microscopic transitions involving the flip¬ 
ping of individual residues between helical and extended conformations. Below we focus on the 
Zimm-Bragg model, as this has been used to characterize the experimental data"^^’^°. In particu¬ 
lar, the elongation parameter s corresponds to the equilibrium constant for adding a single helical 
hydrogen bond. 

We fitted a thermodynamic model to the elongation free energy, here defined as AGei{P, T) = 
—RT In s(P, T). The model includes changes of enthalpy AHq, entropy AS'o and reaction volume 
A Vo under standard conditions, as well as constant differences in heat capacity ACp, isothermal 
compressibility A/3 and linear thermal expansion coefficient Aa to describe the temperature- and 
pressure-dependence. The same model was fitted to the data for RT In s{P,T) reported by Ima- 
mura et al, based on FTIR measurements"^^. The fits to the raw data are shown in Fig. S2^’. In 
Fig. 2 we show the stability diagram for each force field and for experiment; the fitted parameters 
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are listed in Table 1. 


Overall, the stability diagram obtained with the Amber ff03w force field is very similar to ex¬ 
periment, bearing in mind that the experiments were done on a different alanine-based peptide 
(AK20) - we note that the stability diagrams for the AK16 peptide obtained by Hatch et with 
Amber ff03* and TIP3P are qualitatively very similar to those we obtain for Ac-(AAQAA) 3 -NH 2 
with the same force field and water model. This is particularly true in the range of temperature 
and pressure probed by the experiments, indicated by the broken lines in Fig. 2. In both cases, 
increasing pressure clearly stabilizes helical states, as also reflected in the negative reaction vol¬ 
ume of —0.8 cm^.mol”^ for experiment and —1.5 cm^.mol”^ for Amber ff03w (Table 1). It also 
captures quite well the overall enthalpy and entropy changes for adding a helical residue, as previ¬ 
ously noted^^. In contrast, the stability diagram for the ff03* force field differs in some important 
respects. Application of low pressures will have little effect on, or slightly destabilize helical, 
reflected in the small positive AVq = 0.4 cm^.mol”^ (Table 1) for helix formation. Additionally, 
the changes in enthalpy and entropy are almost half the experimental estimates, indicative of too 
low a cooperativity of helix melting^"^. For experiment and both force fields, the changes in heat 
capacity are small. This is in contrast to protein folding^^, and may indicate a limited role for the 
hydrophobic effect in stabilizing helices. For both force fields, the difference in thermal linear 
expansion coefficients, Aa, is small and positive, ~ 4 cm^.mol~^.K“^ This implies a change 
of reaction volume with temperature which will slightly increase the positive reaction volume for 
Amber ff03*, but is insufficient to change the sign of the negative reaction volume using Amber 
ff03w. Over the the range of temperature where liquid water is stable the effect is too small to 
qualitatively change the pressure-dependence of the two models. 

The differences in enthalpy and entropy between the two solvent models have been explained 
in terms of the strength of solvent interactions with the peptide chain^^. These differences also 
result in a more expanded unfolded state for ff03w relative to ff03*. How can the differences in 
reaction volume be understood? The first explanation would be in terms of changes in the peptide 
free energy surface - i.e. different conformations are preferred by each force field, particularly 
for non-helical states. This is undoubtedly a contribution, since these free energy surfaces are 
evidently different. In Fig. 3, we show representative two-dimensional free energy surfaces as a 
function of the radius of gyration and the number of helical residues. As anticipated, the radius 
of gyration of the helical states is similar, but the non-helical states are much more collapsed in 
the case of ff03* compared to ff03w. If one supposes that a more collapsed unfolded state is 
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associated with a smaller volume, then this would be consistent with the observed differences in 
reaction volume for helix formation for the two force fields. 

In order to investigate the above hypothesis, we determined the approximate volume excluded 
to water by each peptide conformation by using the Connolly volume: this is the volume which 
is inaccessible to a sphere of radius 0.14 nm. This is the simplest way in which the volume of a 
configuration can be estimated. The average Connolly volumes for configurations with the same 
number of helical residues are shown in Fig. 4, for the replica at 298 K. The results clearly show 
that more helical states have a smaller total volume for both ff03* and ff03w, which is qualita¬ 
tively consistent with the negative change of volume for helix formation observed experimentally. 
However, it does not explain why under low pressure conditions the reaction volume for TIP3P 
may be positive, and the overall changes of 10-15 cm^.mol”^ per helical residue are about an order 
of magnitude larger than the volume changes estimated from the thermodynamic fits in Table I. It 
also does not explain the reduction of reaction volume for ff03* with pressure, such that it becomes 
negative at sufficiently high pressures. Although we find that helical states do have smaller Con¬ 
nolly volume above 4 kbar, in agreement with earlier findings based on the radius of gyration^^, we 
also find an even greater reduction in volume for non-helical states so that the difference in volume 
between non-helical and helical states would become more positive with increasing pressure, the 
opposite trend to that observed for the total reaction volume of the system. Therefore, a simple 
picture based only on properties of the peptide configurations does not tell the whole story. 

Naturally, the surrounding solvent may also play a role in determining the dependence of sys¬ 
tem volume on peptide conformation, and previous studies have suggested that peptide solvation 
changes with increasing pressure^^’^^, and changes in water structure with increasing pressure are 
known to alter the hydrophobic effecC"^. We probed for this by randomly drawing configurations 
from the 298 K replica from the 1 bar Amber ff03* REMD simulation, and then determining 
for each of these the average system volume using different force fields and system pressures (1 
bar and 4 bar). That is, we effectively remove the influence of different free energy surfaces with 
different force fields by always considering the same set of peptide configurations. By running suf¬ 
ficiently long constant pressure simulations with each of these configurations, restrained to their 
initial position, we can accurately determine the average system volume, as a function of the num¬ 
ber of helical residues. The results of these simulations for the Amber ff03* and ff03w force fields 
are summarized in Fig. 5. We find that this simple analysis captures the reaction volume effects 
inferred from the thermodynamic fits. Namely, at 1 bar pressure, the reaction volume for helix 
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formation is positive with ff03* and negative with ff03w (Fig. 5A,C), while at 4 bar pressure the 
ehange of volume upon helix formation beeomes negative even for TIP3P (Fig. 5B,D). The large 
seatter in the individual system volumes (black data points) indicates that other factors besides 
helicity are important for determining system volume. Nonetheless, when the volume of a given 
configuration with ff03w is subtracted from that with ff03*, these effects are largely eliminated, 
leaving helicity as the main determinant of the difference between the two water models (Figure 
5E-F). This result highlights the importance of solvation in determining the difference in apparent 
volume between the helical and less helical states. 

Although Amber ff03* with the TIP3P water model may fail to describe the correct qualitative 
pressure-dependence of helix stability, it should be noted that the differences in reaction volume 
under consideration are extremely tiny: a change of volume AV = 1 cm^.moP^ per residue, or 
1.66 X 10“^ nm^ per molecule per residue can be related to a change of apparent helix radius AR 
using a helix pitch of AL ~ 0.15 nm and helix radius R ~ 0.45 nm via AR ^ AV/2nRAL, 
yielding a change of apparent radius AR ^ 0.0039 nm. While this is a simplified model calcula¬ 
tion, it serves to illustrate the subtlety of the effect that must be captured. 


IV. CONCLUSION 

The dependence of globular protein stability on pressure has been found to be mainly deter¬ 
mined by cavities in the folded structure, masking other possible effects of interest. In the case 
of helix formation, in fact the opposite trend is found to that for protein folding, namely pres¬ 
sure stabilization of helices. We have shown here that simulations with an accurate water model, 
TIP4P/2005, are capable of capturing the pressure dependence of helix formation. Further, in 
agreement with earlier work, we find that using the TIP3P water model leads to pressure-induced 
destabilization of helices, a qualitatively incorrect result. We further show here that this difference 
is not merely due to the different peptide conformations sampled with that water model. Instead 
there is a critical role for solvent structure in determining the reaction volume changes. Taken 
together, our results emphasize the importance of using an accurate water model for capturing the 
folding of peptides and proteins under different thermodynamic conditions. 

The experimental results obtained for helix formation have proved to be a stringent test for 
simulation models. In future, it would be very interesting to compare simulation results with ex¬ 
perimental data for other model peptides such as /3-hairpins, should such data become available. 
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Experimental kineties results for the pressure-dependenee of helix^° and protein folding®^’^° kinet- 
ies are also a rieh souree of information for future detailed eomparison with moleeular simulation. 
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V. TABLES 


TABLE I. Parameters for fits of helix elongation free energy AGei to thermodynamie model (Eq. 2). 
Numbers in braekets are the error in the last signifieant figure esfimafed by boofsfrap Monfe Carlo. 


Paramefer 

Unifs 

Experimenf 

ff03w 

ff03* 

AHo 

kJ.mol”^ 

-5.2(3) 

-4.1(1) 

-2.9(1) 

ASo 

J.mol“^.K“^ 

-15(1) 

-13.4(2) 

-9.7(3) 

ACp 

J.mol“^.K“^ 

0(1) 

2(1) 

-1(1) 

AEo 

em^.mol”^ 

-0.8(1) 

-1.5(1) 

0.4(2) 

A/3 

xlO“^ em^.mol“^.bar“^ 

6(3) 

22(2) 

-13(2) 

Aa 

xlO“^ em^.mol“^.K“^ 

-3(2) 

4.1(6) 

3.4(4) 
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VI. FIGURES 



Temperature / K 


Bv 


_X 0.8 h 

X 


c 

o 


o 

ns 


■ 

1 bar 


2kbar 

▲ 

4 kbar 

▼ 

8 kbar 

♦ 

12 kbar 



I I I I I I I I I I I I I 

280 300 320 340 360 380 400 420 440 460 480 500 




Temperature / K 


FIG. 1. Pressure-dependence of fraction helix. (A) Amber ff03w, TIP4P/2005 water; (B) Amber ff03*, 
TIP3P water. 


17 
















^ . . 

250 300 350 400 450 


Temperature / K 

FIG. 2. Stability diagrams for helix elongation. The free energy associated with elongating a helix by 
one residue, AGei is obtained from the Zimm-Bragg®^ elongation parameter s as AGgi = —RTlns, as 
a function of temperature and pressure. Shown are the stability diagrams obtained by fitting Eq. 2 to (A) 
experimental RT In s{P,T) obtained from Tmamiira and Katu^^ for the peptide Ac-AA(AAKAA) 3 AAY- 
NH 2 ; (B) constant-pressure replica exchange simulations of the peptide Ac-(AAQAA) 3 -NH 2 with the ff03w 
force field^® and TIP4P/2005 water model^^; (B) constant-pressure replica exchange simulations of Ac- 
(AAQAA) 3 -NH 2 with the ff03* force field^^ and TIP3P water model^^. Dashed magenta box indicates 
approximate region covered by experimental data. Energy units are kJ.rnol”^. 
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FIG. 3. Free energy surfaces for radius of gyration and number of helical residues. (A) Amber ff03w, (B) 
Amber ff03*. Magenta line indicates mean radius of gyration for a given number of helical residues. 
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FIG. 4. Dependence of Connolly volume on helicity. The Connolly volume averaged over configurations 
with the same number of helical residues is shown for (A) Amber ff03* and (B) Amber ff03w. 
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FIG. 5. System volumes determined for a common set of peptide configurations. (A) Amber ff03*, 1 bar; 
(B) Amber ff03*, 4 kbar; (C) Amber ff03w, 1 bar; (D) Amber ff03w, 4 kbar; (E), (F) difference between 
volumes using Amber ff03w and Amber ff03* at 1 bar and 4 kbar respectively. Black data points are 
the average volumes determined for individual configurations (20 per number of helical residues) and red 
symbols are average system volumes for all configurations with the same number of helical residues. 
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FIG. S1. Alternative definition of fraction helix. The temperature-dependent helix formation from REMD at 
constant pressure is shown for (A) Amber ff03w and (B) Amber ff03*. Solid symbols are for helix fraction 
determined from torsion angles (see main text) and open symbols with matching colors are for helix fraction 
determined from DSSP’ using /helix = ?^a/^res^ where is the number of helical residues from DSSP 
and rires = 15 is the number of residues in the peptide. All other details are as in Fig. 1 in the main text. 
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FIG. S2. Fit of thermodynamic model to raw data. (Top) Experimental data from Imamura and Kato'; 
(Center) Simulation data for Amber ff03w; (Lower) Simulation data for Amber ff03*. Symbols are data, 
lines with corresponding colors are fits to the thermodynamic model. 
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